Compensation of decoherence from telegraph noise by means of bang-bang control 
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With the growing efforts in isolating solid-state qubits from external decoherence sources, the 
origins of noise inherent to the material start to play a relevant role. One representative example 
are charged impurities in the device material or substrate, which typically produce telegraph noise 
and can hence be modelled as bistable fluctuators. In order to demonstrate the possibility of 
the active suppression of the disturbance from a single fluctuator, we theoretically implement an 
elementary bang-bang control protocol. We numerically simulate the random walk of the qubit 
state on the Bloch sphere with and without bang-bang compensation by means of the stochastic 
Schrodinger equation and compare it with an analytical saddle point solution of the corresponding 
Langevin equation in the long-time limit. We find that the deviation with respect to the noiseless 
case is significantly reduced when bang-bang pulses are applied, being scaled down approximately 
by the ratio of the bang-bang period to the typical flipping time of the bistable fluctuation. Our 
analysis gives not only the effect of bang-bang control on the variance of these deviations, but also 
their entire distribution. As a result, we expect that bang-bang control works as a high-pass filter 
on the spectrum of noise sources. This indicates how the influence of l//-noise ubiquitous to the 
solid state world can be reduced. 

PACS numbers: 03.65.Yz, 03.67.Lx, 05.40.-a 



In order to implement solid-state quantum information 
processing devices, the decoherence acting on the quan- 
tum states has to be carefully understood, controlled and 
eliminated. So far, research has concentrated on decou- 
pling from external noise sources (like thermal heat baths 
and electromagnetic noise). With the success of this ef- 
fort, noise sources intrinsic to the material such as defect 
states increase in importance and have to be controlled 
in order to improve coherence even further. 

Most external noise sources are composed of extended 
modes in the thermodynamic limit close to equilibrium 
such that their fluctuations are purely Gaussian. Thus, 
their influence can be modelled by an oscillator bath, see 
e.g. 0. However, there are physical situations when this 
assumption fails In particular, this is true for lo- 

calized noise sources with bounded spectra as they occur 
in disordered systems for hopping defect states . Phys- 
ical examples for this situation are background charges 
in charge qubits 0, IE or traps in the oxide layers 
of Joscphson tunnel junctions ||. Such localized noise 
sources are more realistically represented by a collection 
of bistable fluctuators 0,0 (henceforth abbreviated bfls), 
as their noise spectrum is considerably non-Gaussian. If 
many of these noise sources with different flipping times 
are appro priately superimposed, they lead to 1// noise 
IlllOllUl- With the progress of fabrication technology 
and miniaturization of qubits, we expect however that 
there might only be a few fluctuators in a qubit ||. 

We analyze the impact of a single fluctuator in the 



semiclassical limit, where it acts as a source of telegraph 
noise. We apply an open loop quantum control tech- 
nique, namely quantum bang-bang |l2t Hill 1 1 ll] , which is 
designed suitably for slowly fluctuating noise sources. We 
simulate the noise-influenced qubit dynamics with and 
without bang-bang correction by integrating the time- 
dependent Schrodinger equation for each specific realiza- 
tion of the noise. We present the resulting random walks 
around the unperturbed signal on the Bloch sphere and 
analyze the quality of this suppression by an compari- 
son of the ensemble averaged deviations of these random 
walks with and without bang-bang correction. 

We describe our system by the effective Hamiltonian 
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(t) = haZ bfi (t)a% (2) 



where a denotes the coupling strength between the fluc- 
tuator and the qubit and £bfl(i) represents a symmet- 
ric telegraph process that is flipping between ±1, whose 
switching events are Poisson distributed with a mean sep- 
aration Tbfl between two flips. 

On a microscopic level, such noise is typically gener- 
ated by coupling the qubit to a two-state impurity, which 
is in turn coupled to a heat bath causing the two-state 
system to flip randomly and incoherently. Our model 
corresponds to the semiclassical limit and should be ac- 
curate whenever the coupling of the impurity to the bath 
is much stronger than its coupling to the qubit 0, such 
that the qubit does not act back on the noise source. 
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The assumption of a symmetrical telegraph process cor- 
responds to a high bath temperature compared to the 
impurity level spacing. This restriction is not essential 
for the following investigations, for an asymmetric noise 
signal would only produce an additional constant drift. 

We describe the resulting evolution of the noise influ- 
enced qubit by a stochastic Schrodinger equation 0, 
with the time-dependent Hamiltonian (J2J. For any initial 
state of the qubit, we numerically integrate 



tp(t) = Texp( -i/h \ Hf(t')dt')il)(0) (3) 



viation from the unperturbed signal 



A ^w = J^Efe(^^»w) 2 ( 4 ) 



In order to characterize the degree of noise suppression 
by means of bang-bang control, we define the suppression 
factors for a given time tg 
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The deviation as a function of time is plotted in Fig. Q 



with T the time-ordering operator and ip the state vector 
of the two-state system. The result is a random walk 
on the Bloch sphere, which is centered around the free 
precession corresponding to e q and A q . 

We implement the following idealized open quantum 
control scheme: apply an infinite train of 7r-pulses on the 
qubit with negligibly short pulse durations and a con- 
stant separation time Tbb between neighboring pulses. 
In doing so, we intend to average out the a z parts of 
the effective Hamiltonian (and thereby in particular the 
noise term) on time scales large compared to Tbb- This 
is accomplished by iteratively spin-flipping the qubit and 
thus effectively switching the sign of the noisy part of the 
Hamiltonian. This mechanism thus works analogously 
to the well-known spin-echo procedure, specifically the 
Carr-Purcell procedure of NMR We expect to com- 
pensate a fraction of the telegraph noise effects: the size 
of the random walk induced by the noise is determined 
by the typical time separation of the fluctuators influ- 
ence between two flips Tbfl and its coupling strength a 
and scales roughly with ar^a [jj. Using bang-bang, the 
bfls influence remains uncompensated for at most a sin- 
gle bang-bang period. Thus, we reduce the influence of 
the bfl randomly by an average factor of ibfl/rbb- 

As generic conditions of the system dynamics we con- 
sider for the numerical simulations e q = A q = 51o- With- 
out loss of generality, we assume (a%) = +1 as an initial 
state. If there were no noise, the spin would precess on 
the Bloch sphere around the rotation axis a% + a%. So we 
expect for not too large an interaction strength (a <C 1) 
a slight deviation of the individual quantum trajectory 
from the free evolution case. We take a = 0.1 for our 
coupling strength. All the following time and energy 
measures are given in units of the unperturbed system 
Hamiltonian: our time unit is rg ys — 1 / S7o and our en- 
ergy unit is AE = + A^ = v^o- Note that in 

these units, a period lasts 7rr sys /-\/2. We have integrated 
Equ. © and averaged over N ~ 1000 realizations. The 
time scale ratio Tbfl/Tbb = 10 if not denoted otherwise. 
We characterize our results by the root-mean-square de- 




FIG. 1: Time evolution of the mean deviations for bfl-induced 
random walks with and without bang-bang. The straight lines 
are square-root fits of the analytical derived random walk 
model variances (plotted as triangles). Inset: Transverse and 
perpendicular components of bang-bang suppressed noise. 



We recognize that the total deviations on intermediate 
time scales are suppressed by a ratio of ~ 10. Detailed 
analysis shows that the tangential and the orthogonal 
deviation, corresponding respectively to dephasing and 
relaxation, are of the same size for the uncompensated 
case. In contrast, the bang-bang modulation mostly com- 
pensates the dephasing-type deviation, as shown in the 
inset of Fig. {IJ. 

We now develop analytical random walk models for 
our system. Although the random walk on the Bloch 
sphere is in general two-dimensional, bang-bang control 
effectively reduces it to a one-dimensional model, rep- 
resenting the relevant perpendicular part. We restrict 
ourselves to the long-time limit. 

For simplicity, we replace the fluctuating number of 
random walk steps for a given time At of noisy evolu- 
tion by its expectation value At/r^n [ijjj. This allows to 
use the number of random walk steps as time parameter. 
We encounter different one-step-distributions, depending 
on whether the number of steps is odd or even, corre- 
sponding to an "up" or "down" state of the bfl |22| . The 
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step-size distribution of the bfl model in our small devi- 
ation regime is given from Poisson statistics 



ff,bfl ( \ 

^odd/cvcnV J V 







(6) 



with = arbfl as a typical random walk one-step devia- 
tion. r un it is a time unit, corresponding to a discrete step 
length iE un it = a r un it of the random walk. 6(x) denotes 
the Heaviside step function. We neglect the correlations 
between transverse and perpendicular deviations as they 
average out in the long-time limit. 

For the bang-bang suppressed random walk, the flip- 
ping positions of the bfl-noise sign in the bang-bang 
time-slots are essentially randomly distributed as long as 
7"bb <^ T bfl- We find a homogenous step-size distribution 
between zero deviation and a maximum 7 
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The ubiquitous occurs, because the bang-bang se- 
quence also averages over the static e q -term and hence 
slows down the free evolution. 

By means of these one-step probability distributions, 
we are able to calculate via convolution the distributions 
for 27V-step random walks. Specifically, they are the in- 
verse Fourier transforms of the TV-fold products of the 
Fourier transforms of the two-step distribution ff^]. For 
the uncompensated case, we find 
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whereas for the compensated case 
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Already for random walk step-numbers in the order 
of ten, the resulting distributions are almost Gaussian. 
Their standard deviations give the rms deviations of the 
random walk models plotted in Fig. QJ. As expected, 
they grow as a square-root of the number of steps. 

The above integrals can be evaluated analytically using 
the saddle point approximation. We find variances of 



CTbfi(iV) = V2AW bfl 
for the pure bfl random walk and 



(10) 



(11) 



for the compensated one. In the large- iV limit, this model 
shows excellent agreement with the simulation. 

Beyond predicting the variance, our analysis also al- 
lows evaluation of the full distribution. We compared 



evolution with and without bang-bang compensation via 
simulations with 10 4 realizations and calculated the full 
distribution function for a evolution time to = Ts ys . The 
numerical histograms of the deviation with their respec- 
tive one- and two-dimensional Gaussian fits are shown in 
Fig.©. 
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FIG. 2: Histograms of the deviation from free evolution with 
and without bang-bang and fits to the expected two- (pure 
bfl case), respectively one-dimensional (bang-bang corrected 
evolution) random walk statistics. Numerical data collect 
over 10000 realizations at a fixed time to = rs ys . With 
Tbfl = 0.01rsy S the random walk distributions are calculated 
for N = rsys/Tbfl = 100 steps. (NB: The :r-axis scale of the 
right graph is 15 times smaller than that of the left graph.) 



We observe that not only the bang-bang compensated 
distribution is much narrower than the uncompensated 
distribution, but also that its shape is qualitatively dif- 
ferent: its maximum is at zero error whereas the uncom- 
pensated distribution has its maximum at a finite error 
I Act I « 0.01 and zero probability of zero error. 

We have systematically studied suppression factors for 
different ratios of the switching time t^/t^ at a con- 
stant fluctuator flipping rate Tbfl = 10 _2 T sys and evolu- 
tion time to = r sys- The numerical data in Fig. (3) show 
that the suppression efficiency is linear in the bang-bang 
repetition rate, S = /-tTba/rbb- The numerically derived 
value of the coefficient, /^numerical ~ 1-679, is in excellent 
agreement with our analytical result /J. a naiyticai = v3 ~ 
1.732 from the saddle point approximation, Equs. ifTTijl 
and l|ll|) . This small discrepancy reflects the correlations 
between the transverse and longitudinal random walk in 
the uncompensated case, see Fig. 

We have demonstrated the ability of a bang-bang pro- 
tocol to compensate environmental fluctuations with fre- 
quency u> <C 1/Vbb- Thus, bang-bang is acting as a "high 
pass filter" for noise with a roll-off frequency of 1/Vbb- 
Evidently, the bang-bang correction is suitable for sup- 
pressing the impact of telegraph noise on qubits and can 
enhance the coherence by orders of magnitude. The ap- 
plication of the scheme which we outlined requires a rel- 
atively strict separation of time scales: One has to be 
able to flip the spin very rapidly, typically two orders 
of magnitude faster than Tbfl. It remains to be investi- 
gated how this scheme works with pulse durations that 
are finite rather than infinitesimal. Moreover, we have 
assumed that the environment produces symmetric tele- 
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FIG. 3: The suppression factor <S t() (r b fl/r b b) = ^bb 8 ^ ) eva l" 
uated for to = rs ys as a function of the ratio of the nipping 
time TbH and the bang-bang pulse separation Tbt,- 

graph noise regardless of the qubit dynamics. Clearly, 
the issue of when one may neglect feedback effects be- 
tween the qubit and bfl must be critically revisited in the 
low-temperature limit. We speculate that the setup is 
promising for l//-noise, as in particular the most harm- 
ful and predominantly low-frequency fraction of a corre- 
sponding ensemble of fluctuators would be compensated 
most strongly. Finally one has to be aware that also 
the static term of the Hamiltonian is averaged out, and 
this generally reduces the degree of control on the qubit. 
This is only a technical constraint, however, as one could 
imagine interchanging two different types of bang-bang 
pulses to admit corresponding quantum gate operations. 

Another approach for decoupling from slow noise is 
to choose an appropriate working point with a domi- 
nant term £la x in the static Hamiltonian. The action 
of this term can be understood as a rapid flipping of the 
spin, similar to our bang-bang protocol. Using a Gaus- 
sian approximation the noise from the bfl with standard 
rate expressions (e.g., 0), it can be shown that the de- 
phasing rate reads = ev/rbflfi 2 instead of — evT b fl, 
which corresponds to the same amount of reduction as 
in our case. This scheme has been implemented in su- 
perconducting qubits [13] ■ In that case, it turned out 
that because the <j x term was limited by fabrication, this 
consideration led to a major redesign. Our compensation 
scheme purely relies on external control and thus keeps 
the hardware design flexible. 

A related problem has been addressed in Ref. |2lj . 
which deals with bang-bang suppression of Gaussian 1//- 
noise, i.e., a bosonic bath with an appropriate sub-Ohmic 
spectrum. That system is treated in the weak-coupling 
approximation, i.e. it assumes S(u>)/u/ < 1 at low fre- 
quencies where S(ui) is the noise spectral function. Both 
assumptions are serious constraints in the 1 / /-case 0, |^ . 
Our work is not constrained to weak coupling, takes the 
full non-Gaussian statistics of telegraph noise into ac- 
count, and gives the full resulting distribution of errors. 

In summary, we examined the decoherence of a sin- 



gle qubit from a single symmetric telegraph noise source 
and proposed an adequate open quantum control com- 
pensation protocol for suppressing its impact. We simu- 
lated the qubits dynamics using a stochastic Schrodinger 
equation and analyzed its deviation from free evolu- 
tion. We formulated analytically solvable one- and two- 
dimensional random walk models, which are in excellent 
agreement with the simulations in the long time limit. 
Specifically, we show quantitatively, how the degree of 
noise compensation is controlled by the ratio between bfl 
flipping time scale and bang-bang pulse length. We give 
the full statistics of deviations in both cases. 
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